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Abstract 

We report on a first-principles study of the structure of O/Rh(001) at half a 
monolayer of oxygen coverage, performed using density-functional theory. We 
find that oxygen atoms sit at the center of the black squares of a chess-board, 
c(2 x 2), pattern. This structure is unstable against a rhomboid distortion of 
the black squares, which shortens the distance between an O atom and two of 
the four neighboring Rh atoms, while lengthening the distance with respect to 
the other two. We actually find that the surface energy is further lowered by 
allowing the O atom to get off the short diagonal of the rhombus so formed. 
We predict that the latter distortion is associated with an order-disorder tran- 
sition, occurring below room temperature. The above rhomboid distortion of 
the square lattice may be seen as a rotation of the empty, white, squares. Our 
findings are at variance with recent claims based on STM images, according to 
which it is instead the black squares which would rotate. We argue that these 
images are indeed compatible with our predicted reconstruction pattern. 
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I. INTRODUCTION 



The study of the physical and chemical properties of the (001) surface of rhodium is 
an important step towards the ambitious goal of understanding the catalytic properties of 
transition metals. Recently, the clean Rh(001) surface has attracted experimental [l] |3[ 
and theoretical |4]-|9| interest because of the anomalously small interlayer relaxation at the 
surface, an unexpected result which is in contrast with the general behavior of transition- 
metal surfaces. One more reason of interest is due to the behavior of the surface upon oxygen 
adsorption. Oxygen adsorption on Rh (001) is known to be dissociative and to saturate at 
half a monolayer, independently on the adsorption temperature. At this coverage a 2 x 2 
reconstruction has been observed by SPA-LEED and confirmed by STM [[l(|. The oxygen 
atoms sit in the troughs formed by four first-layer rhodium atoms and fill these sites in 
a c(2 x 2) geometry. This geometry may be seen as a chess-board whose black squares 
are occupied by oxygen atoms, while the white ones are empty. Within this picture, the 
reconstruction observed in Ref. ]TU| has been described as a rotation of the black squares, 
resulting in a (2 x 2)p4g symmetry (see Fig. |I|a). This distortion preserves the shape of the 
black squares, while the white ones become rhomboid. A similar behavior is observed for 
nitrogen and carbon adsorbed on the (001) surface of nickel, where the rotation angle of the 
squares is much larger and the clock reconstruction more evident |TT|-|nj . 

In this paper we present an ab-initio study of the oxygen adsorption on the Rh(001) 
surface. In agreement with experimental findings, we find that the preferred adsorption site 
is the hollow one, at the center of the square formed by four first-layer rhodium atoms. 
However, at variance with recent claims based on STM experiments |Kj, we find that it is 
the white squares which rotate rather than the black ones (see Fig. |l]b). In spite of this ap- 
parent discrepancy, we argue that our predicted reconstruction pattern is indeed compatible 
with experimental findings. Actually, the symmetry of our predicted reconstruction pattern 
slightly differs from what appears in the experimental STM pictures, in that ad-atoms get 
off the center of the distorted squares, thus resulting in an asymmetric clock reconstruction. 
The existence of two inequivalent low-symmetry sites that oxygen atoms can occupy in the 
rhombi should give rise to an order-disorder transition whose properties, however, lie beyond 
the scope of the present study. The amplitude of the asymmetric component of the recon- 
struction might be too small to be detected experimentally. Furthermore, it is likely that the 
temperatures at which these STM images have been taken lie above the critical temperature 
of transition. In order to understand the driving mechanism of the reconstruction, this is 
better visualized as a rhomboid deformation of the squares filled by adsorbed oxygen atoms. 
This distortion can then be rationalized by a simple chemical model of the re-bonding of O 
atoms at the surface. 



II. COMPUTATIONAL METHOD 



Our calculations are based on density functional theory within the local-density ap- 

|18[. The 



proximation []16 , |17 |, using Ceperley- Alder exchange-correlation (XC) energies 
one-particle Kohn-Sham equations are solved self-consistently using plane-wave basis sets 
in a pseudopotential scheme. Because of the well known hardness of the norm- conserving 
pseudopotentials for the O and — to a lesser extent — Rh atoms, we make use of ultra-soft 
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pseudopotentials |1| which allow an accurate description of the and Rh valence pseudo- 
wave-functions with a modest basis set including plane waves up to a kinetic-energy cutoff 
of 30 Ry. In the case of Rh, we found it convenient to treat the s and p channels using a 
norm-conserving potential, while the ultra-soft scheme is applied only to the hard d orbital 



22]. Brillouin-zone integrations have been performed using the Gaussian-smearing [23 
special-point |24| technique. We find that the structural properties of bulk rhodium are well 
converged using a first-order Gaussian smearing function [f23[ of width a = 0.03 Ry and 10 
special k-points in the irreducible wedge of the Brillouin zone (IBZ). The isolated surface is 
modeled by a periodically repeated super-cell. We have used the same super-cell for both 
the clean and the O-covered surfaces. For the clean surface we have used 5 atomic layers plus 
a vacuum region corresponding to ~ 6 layers. For the O-covered surface, the 5 Rh layers 
are completed by one layer of O molecules on each side of the slab: in this case the vacuum 
region is correspondingly reduced to ~ 5 atomic layers. We have used the same Gaussian- 
smearing function as in the bulk calculations with a (12, 12, 2) Monkhorst-Pack mesh P^j 
resulting in 21 special k-points in the lxl surface IBZ. Convergence tests performed with 
a value of a twice as small and a correspondingly finer mesh of special points resulted in 
no significant changes in total energies and equilibrium geometries. The latter are found by 
allowing all the atoms in the slab to relax until the force acting on each of them is smaller 
than 0.1 x 10~ 3 Ry/a . 



III. RESULTS 

For the clean surface we find that — in agreement with other ab-initio calculations [f§-§J- 
the first layer relaxes inward by 3.8 ±0.2 % (see Table |), while the second layer is practically 
unrelaxed (the error is essentially a finite-size effect, and it is estimated by repeating the cal- 
culation using thicker slabs). As it was noted in the introduction, the value of the first-layer 
relaxation reported in the experimental literature is anomalously small (the first interlayer 
spacing is practically equal to its bulk value within error bars). Recently, Cho and Scheffier 
H pointed out that a proper account of the vibrational contribution to the surface free 
energy may results in a reduction of the inward relaxation of the first layer, thus bringing 
theoretical predictions in better agreement with experiments. No attempts of estimating 
these vibrational effects have been done in the present work. Our calculated values for the 
surface energy, a, and work function, <f), are a ~ 1.4 eV/atom and (f) ~ 5.5 eV (see Table |). 

A special care must be payed when estimating the latter, for the XC potential goes to 
zero rather slowly in the vacuum region. As the electrostatic (Hartree) potential converges 
much more rapidly than the XC one to its vacuum value, it is convenient to evaluate the 
vacuum level by simply neglecting the XC contributions. As a preliminary study of oxygen 
adsorption, we performed a couple of simple ab-initio molecular dynamics simulations of an 
O2 molecule impinging onto the surface. We have tried two possible initial conditions for 
the molecule, in both of which the axis of the molecule is parallel to the surface. In the 
first case the projection of the center of the molecule on the surface falls on top of a surface 
atom, while in the other it falls on a bridge site. The initial velocity of the molecule is 
orthogonal to the surface and its modulus corresponds to a temperature of 300°K, while the 
surface is initially assumed to be at zero temperature. The initial distance of the molecule 
from the surface is ~ 6 a.u. which is essentially in the vacuum. When the molecule arrives 
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at the surface the molecular bond breaks and the surface heats up. In order to find the 
ground-state atomic configuration, we have slowly cooled down the system, ending in the 
structure displayed in Fig. |2|: the oxygen atoms adsorb in every second hollow site of the 
surface, resulting in a c(2 x 2) structure, where the squares whose centers are occupied by 
an ad-atom are deformed into rhombi; the oxygens do not stay in the middle of the short 
diagonal of the rhombus so formed, but move away in the orthogonal direction, becoming 
essentially three-fold coordinated. As a result of the reduced symmetry so obtained, the 
atomic rows of the first Rh layer are no longer equivalent to each other: every second row 
is formed by atoms which have 2 oxygen neighbors, while the others have only 1 neighbor. 
As a consequence the first rhodium layer results to be buckled, the 2-fold coordinated rows 
leaning « 0.08 A outwards. 

This asymmetric clock reconstruction can be imagined as the succession of two steps: in 
the first, the occupied, black, squares distort while the ad-atoms stay at the center of the 
rhombi (symmetric clock reconstruction) ; in the second, the ad-atom depart from the centers 
of the rhombi, in order to achieve three-fold coordination. Both the unreconstructed and 
symmetric clock structures result to be saddle points whose energies are reported in Table 
|. We have explicitly verified that the unreconstructed structure is stable against analogous 
distortions affecting the white squares. For the sake of completeness, in Table | we also 
report the relative surface energies corresponding to two other adsorption sites (see Fig. 
|3|), which are however much larger than the energy differences between the unreconstructed 
hollow structure and the corresponding clock reconstructions. 

The work function is larger at coverage = ~ than for the clean surface (see Table 
|), thus indicating that electrons tend to transfer from the Rh substrate to the oxygen 
layer, so as to increase the surface dipole. In Fig. | we display the adsorption energy, E ad , 
as a function of coverage, 6, defined as: E ad {Q) = er(O) — er(0) — ®£^o 2) where cr(6) is 
the surface energy of the oxygen-covered system, and Eq 2 is the energy of the isolated 2 
molecule. Coverage O = 1 is realized by filling all the hollow sites; for G = | and = |, 3 
and 1 of the 4 possible hollow sites in a 2 x 2 substrate super-cell have been occupied. The 
minimum energy corresponds to a coverage in the range \ < < 1. In this ranee, however, 
the curve is very flat and the variation of the adsorption energy (« 100 meV) is smaller 
than the expected accuracy of the terms which enter its definition. All we can safely predict 
is that the stable coverage state lies somewhere in between = | and = 1. The value 
of the most stable coverage is determined by a trade-off between the adsorption energy of 
an isolated oxygen molecule, which tends to favor a high coverage, and the oxygen-oxygen 
electrostatic repulsion due to charge transfer, which becomes more effective when the average 
— distance becomes smaller than some typical screening length, and which tends instead 
to favor a low coverage. 



IV. DISCUSSION 

From Fig. it is evident that the oxygen sub-lattice forms a zig-zag arrangement which 



is not observed in the experiments | 10|| . However, there are two equivalent threefold sites 
for each cell occupied by an oxygen ad-atom, one on each side of the bridge. Neglecting the 
interactions between different adsorption sites, each of them is therefore two-fold degenerate. 
We postulate that the system undergoes an order-disorder transition at some critical tem- 
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perature which depends on the magnitude of the adsorbate-adsorbate interactions, and that 
this temperature is lower than that at which the STM images have been taken. The fact 
that no disorder in the O ad-layer appears from these images indicates that the frequency 
of oxygen barrier crossing is much larger than the scanning frequency of the STM. This is 
compatible with our estimate of the barrier height (« 30 meV, as obtained by comparing 
the energies of the symmetric clock and the asymmetric clock structures) which is of the 
order of the room temperature. 



The STM pictures reported in Ref. |T0[ are such that only Rh or O atoms are visible 
in turn in a same picture. Because of this, it is difficult to judge whether the squares that 
undergo a rhomboid distortion are those filled by an ad-atom, or they are instead empty. 
What appears to be rather evident is that the distorted squares display a depletion with 
respect to the undistorted ones. In analogy with a similar reconstruction which is observed 



to occur in N/Ni(001) and which is much better characterized [fT5 ] , the authors of Ref. |10 



conclude that O atoms are thrust down the first Rh layer in the middle of the undistorted 



squares. We have simulated STM pictures using the Tersoff-Hamann model [27], but we 



have not been able to eliminate the images of the O atoms from any of our simulated 
pictures. In all of them, O atoms appear as protrusions in the middle of the distorted 
squares, in qualitative agreement with the experimental observation that undistorted squares 
are depleted with respect to the distorted ones. 

In Fig. ||we display the surface-projected densities of states (SDOS) of the clean and of 
the oxygen- covered Rh(001) surface, along with their decompositions into various atomic- 
like contributions. The SDOS of the clean surface differs from its bulk counterpart mainly 
because of its more pronounced narrowness due to the lower coordination of surface atoms 
0,21], and it is almost entirely determined by its d-\ike component (see Fig. |5]a). In the bulk 
FCC structure the three d xy , d xz and d yz are degenerate, and so are d 3z 2_ r 2 and d x 2_ y 2. This 
degeneracy is partially lifted at the surface because the surface atoms lose the neighbors in 
the direction orthogonal to the surface. For the (001) surface, the d xz and d yz orbitals are still 
equivalent by symmetry and, hence, degenerate. The position of the oxygen atomic p-level 
is around 3.1 eV below the Fermi energy of the oxygen- covered surface. Inspection of Fig. 
|5]d shows that the Px = Py level gives rise to a bonding and an anti-bonding main peaks, 
respectively below (—5.8 eV) and above (+1 eV) the Fermi level. Upon oxygen adsorption, 
first-layer Rh atoms become locally inequivalent according to whether the neighboring oxy- 
gen atoms are aligned along the x or y directions (see Fig. ^|e). The projected densities of 
states (PDOS) plotted in Fig. [|d refer to those Rh atoms which have O neighbors aligned 
along the y direction (the PDOS of the other Rh atoms can be obtained by simply exchang- 
ing x with y). It is also easy to recognize that the same bonding and anti-bonding features 
occur in the d yz band whose atomic orbitals have lobes oriented towards the adsorbed oxy- 
gen, while the d xz band remains similar to that of the clean surface (see Fig. ||]c) because in 
that direction the surface sites are empty. The Px orbitals make bonds with the d xz orbitals 
of the Rh atoms along the x direction, while the Py orbitals hybridize with the d yz orbitals 
of the Rh atoms along the y direction. We find that the c(2 x 2) structure is not the most 
stable one and that a (2 x 2)p4g reconstruction occurs which shortens and strengthens the 
bonds O— Rh along one direction, while lengthening them in the perpendicular direction, 
as illustrated in Fig |5|e. The amplitude of the distortion depends rather sensitively on the 
lattice parameter: using our calculated lattice parameter, we estimate the distortion to be 
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~ 4%; if one uses a lattice parameter 1% larger, the amplitude of the distortion also is 
increased, reaching a value ~ 6%. The opposite occurs if the lattice parameter is reduced. 
The tendency of the system to strengthen two of the four rhodium bonds at the expenses of 
the other two results in a net lowering of the surface energy of ~ 3 meV/atom, thus stabiliz- 
ing the (2 x 2)p4g structure. In the asymmetric clock more stable structure the amount of 
the distortion is ~ 11% (using our calculated lattice parameter) and the surface energy gain 
is ~ 30 meV/atom. This shows that the reconstruction occurs because the optimal O— Rh 
bond-length is shorter than that realized in the ideal geometry. The chemical contribution 
to the energy lowering which determines the distortion is illustrated in Fig. ^|f which shows 
the differences between the Rh^ and O p PDOS after and before the reconstruction. In both 
cases, we notice a push of states towards lower energies. 

V. CONCLUSIONS 

Our ab-initio results confirm the experimental evidence of the dissociative character of 
the oxygen adsorption on Rh (001) and that the favored adsorption site is the fourfold one. 
A (2 x 2)p4g reconstruction of the surface is also predicted, in agreement with SPA-LEED 
data. At variance with claims based on recent STM work, we find that this reconstruction 
is due to a rhomboid distortion of the squares formed by first-layer Rh atoms which have 
an O atom adsorbed in the middle, and that this structure is unstable with respect to a 
departure of the ad-atoms from the centers of the resulting rhombi. We argue that the 
experimental evidence upon which these claims are founded is rather weak and we suggest 
therefore that further experimental work is needed to fully characterize the reconstruction 
of this surface. Furthermore, we find evidence that an order-disorder transition should 
occur at some temperature below that at which the STM pictures have been taken. Our 
findings are substantiated by a simple chemical model of the mechanisms responsible for the 
reconstruction. 

Our calculations were performed on the SISSA IBM-SP2 and CINECA-INFM Cray- 
T3D/E parallel machines in Trieste and Bologna respectively, using the parallel version 
of the PWSCF code. Access to the Cray machines has been granted within the Iniziativa 
Trasversale Calcolo Parallelo of the INFM. 
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TABLES 



TABLE I. Structural data for the three oxygenated structures investigated (see Fig. ^) and 
for the clean surface. The coverage is = |. do is the bulk lattice spacing, doi is the distance 
between the oxygen atoms and the first rhodium layer, di2 the distance between the first and the 
second layer. For the asymmetric clock reconstruction the two given numbers refer to the two 
inequivalent first layer rhodium atoms (see Fig. ||). 5 is the amplitude of the movement of the 
first-layer rhodium atoms upon distortion (see Fig. ||(b)), <ft is the work function and a is the 
surface energy. 
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FIGURES 



FIG. 1. Two possible clock reconstructions of the O/Rh(001) surface, resulting in a (2 x 2)pAg 
structure. In the black reconstruction (a) the squares with an O atom in the middle rotate, while 
the empty ones distort to rhombi; in the white reconstruction (b), the opposite occurs. 

FIG. 2. Equilibrium structure of O/Rh(001) as obtained by our simulated-annealing procedure. 
The thick line indicates the unit cell. The thin line indicates the rhombus and its shorter diagonal. 
The oxygen atoms are alternatively shifted orthogonally with respect to the shorter diagonal. The 
atoms of the first surface layer (depicted with a brighter tone) are « 0.08 A higher then the others. 

FIG. 3. Sketch of the various O adsorption sites of the Rh(001) surfaces considered in this 
work. Small dark circles: O atoms. Large lighter circles: first-layer Rh atoms. Large darker 
circles: second-layer Rh atoms. 

FIG. 4. Oxygen adsorption energy, E a d (eV), as function of coverage, 0. 

FIG. 5. a) Surface projected density of states for the clean Rh(001) surface; b) Surface pro- 
jected density of states for the the oxygen covered one Rh(001) surface; c) Projection of the SDOS 
onto the d xz = d yz first-layer rhodium orbitals for the clean surface; d) Projection of the SDOS 
onto the p x = p y oxygen orbitals and of d xz and d yz first-layer rhodium orbitals for oxygenated, 
unreconstructed, surface; e) Sketch of the (2 x 2)p4g, clock, reconstruction. The amplitude of the 
deformation is exaggerated for clarity; f) Differences between the O p — and Rh^— projected density 
of states in the reconstructed and unreconstructed structures. 
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